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Abstract The study of the TeV emission from extra- 
galactic sources is hindered by the uncertainties on the 
diffuse Extragalactic Background Light (EBL). The re- 
cent H.E.S.S. results on the blazars 1ES 1101-232 and 
H 2356-309 represent a breakthrough on this issue. Their 
unexpectedly hard spectra allow an upper limit to be 
derived on the EBL in the optical/near-infrared range, 
which is very close to the lower limit given by the re- 
solved galaxy counts. This result seems to exclude a large 
contribution to the EBL from other sources (e.g. Popula- 
tion III stars) and indicates that the intergalactic space 
is more transparent to 7-rays than previously thought. A 
discussion of EBL absorption effects and further obser- 
vational tests with Cherenkov telescopes are presented. 

Keywords 7-rays • AGN • EBL • extragalactic • diffuse 
background • blazar 

PACS 98.54.Cm • 98.62.Ai • 95.85.Kr • 95.85. Jq • 
95.85.Pw 



1 Introduction 

The diffuse Extragalactic Background Light (EBL) con- 
tains unique information about the epochs of formation 
and the history of evolution of galaxies, since it con- 
sists of the sum of the light produced by all extragalac- 
tic sources over cosmic time. Its Spectral Energy Dis- 
tribution (SED, Fig. [1} is characterized by two broad 
bumps^lj, produced by thermal emission contributed di- 
rectly from stars at optical (Opt) and near infrared (NIR) 
wavelengths, and by dust which partly absorbs and re- 
emits the starlight at longer wavelengths (FIR). The 
Opt-NIR range has gained special interest in recent years 
because some direct measurements between 1 and 3 /im 
have suggested the presence of a significant excess above 
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Fig. 1 Spectral Energy Distribution of the diffuse Extra- 
galactic Background Light. Open symbols correspond corre- 
spond to the integrated light from resolved galaxies, and thus 
represent lower limits for the EBL (refs in [5])- From left to 
right: HST data (green and red sq uare s) , Spitzer data (blue 
squares ), I SO data (red triangles 22 ), Spitzer data (blue 
squares 23 ; small filled squares are the corresponding esti- 
mates for the total galaxy emission). The other filled sym- 
bols show the different EBL estimates (compilation from [T] 
I5lf24] h Above 300 /im the CMB starts to dominate. The black 
line shows the double-humped phenomenological curve used 
as reference shape for the EBL SED (see [I][5]). The upper 
axis shows the energy of the 7-ray photon which interact with 
the corresponding EBL photon at the peak of the 77 cross- 
section. 

the contribution of normal galaxies, which has been in- 
triguingly interpreted as the signature of redshifted UV 
light from heavy Population-Ill star formation at z ~ 10 
[2] (though with severe energy-budget problems, see [3]). 
Unfortunately, direct EBL measurements are subject to 
possible large systematic uncertainties due to the diffi- 
culty of an accurate modelling and then subtraction of 
the bright foregrounds, mainly consisting of zodiacal (in- 
terplanetary dust) light pp. 

An alternative approach [4] is provided by the mea- 
surement of very high energy (VHE) 7-rays from extra- 
galactic sources, through the detection and identifica- 
tion of absorption features in the VHE spectra caused by 
the interaction with EBL photons (77 — > e + e~). High- 
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Fig. 2 Effects of changing the EBL normalization (expressed as scaling factor P: left panel, l.Ox and 0.45x) on the 
attenuation factor e~ T , for a fixed redshift (center panel: full and dashed blue lines, respectively). In the center panel, the 
effect of increasing redshift is also shown (red and blue full lines, respectively) for a fixed EBL SED (P1.0). The attenuation 
curves represent directly the shape of an observed VHE spectrum, if the source has an emission spectrum parallel to the 
upper axis. All calculations were performed assuming Ho = 70 km/s/Mpc, O m = 0.3, Qa = 0.7. Right panel (full lines): the 
difference in photon index between emitted and observed spectrum (when fitted with a power-law model between 0.2 and 2 
TeV) is plotted as a function of redshift, for two EBL normalizations. The dashed line shows the milder deformation of the 
original spectrum in the energy range 1-5 TeV, due to the oc A" 1 EBL slope between 2-3 and 10 /im. 



energy peaked blazars (HBL) represent to this respect 
excellent sources of 7-ray beams, thanks to their promi- 
nent TeV emission and high apparent luminosity, which 
makes them detectable from large distances and/or with 
good photon statistics. The "downside" is that they are 
highly variable and characterized by a very wide range 
of possible spectra. Though a lot has been learned about 
them, the present understanding of their radiation pro- 
cesses is not yet complete enough to reliably predict the 
intrinsic TeV spectra, and thus to disentangle absorp- 
tion from intrinsic features. Conversely, the still large 
uncertainties on the EBL prevent a sufficiently accurate 
reconstruction of their intrinsic spectra in order to con- 
strain the emission models. The studies of blazars and 
EBL are therefore tightly coupled, in a classic case of 
"one equation, two variables" problem [HE]. 

So far, only hints and guesses could be made from the 
TeV detected objects (in particular from the more dis- 
tant 1ES 1426+428 [8] and PKS 2155-304 [9EQ]), since 
the reconstructed spectra could be generally explained 
also as blazars' intrinsic features. The surprising hard- 
ness of the H.E.S.S. spectra of 1ES 1101-232 and H 2356- 
309 [5, 6]) represent a breakthrough in this respect, posing 
the problem of conflict between blazar spectra and EBL 
estimates in a much more severe and compelling way. 



2 EBL effects on TeV 7-ray spectra 

To illustrate the effects of absorption, as reference shape 
for the EBL SED we adopted the phenomcnological curve 
used in jHTS] (labelled P1.0 in Figure 1), which is designed 
to be in general agreement with the EBL spectrum ex- 
pected from galaxy emission JTj (modulo normalization). 
Such EBL SED imprints a specific shape onto the orig- 
inal spectrum, as illustrated in Fig. 2 (center panel): 



above ~0.1 up to ~l-2 TeV it causes a strong steep- 
ening, followed by a flattening between 2 and ~7 TeV, 
where the observed spectrum partially recovers the orig- 
inal slope (owing to the decline approximately a A -1 
above few microns, which causes the optical depth to 
become nearly constant with energy[3]). Further on, a 
sharp cut-off occurs above ~8-10 TeV caused by the EBL 
SED rising again above 10 /zm. 

This yields two consequences: a) with the current 
systematic uncertainties of Cherenkov telescopes (CT), 
if the intrinsic spectrum is a power-law the observed 
spectrum can still be well approximated by a power-law 
model in the range 0.2-2 TeV, but of steeper spectral in- 
dex (r o b s > Tint)] b) measuring a spectrum above ^200 
GeV, no cut-off is expected from EBL absorption up 
to ~6-7 TeV: any cut-off measured by CT in such en- 
ergy range is more likely to be intrinsic (unless consid- 
ering an EBL NIR spectrum much flatter than oc A -1 ). 
Note also that in the range 0.2-10 TeV breaks or cut- 
offs in the emitting particle distribution of blazar jets 
are quite common, as shown by the corresponding syn- 
chrotron emission in the Opt-X-ray range. Such intrinsic 
steepening might cancel out the expected flattening due 
to EBL absorption. In other words, the absence of such 
flattening in the TeV spectra of some blazars does not im- 
ply an evidence for a different EBL spectrum. To detect 
such feature, a power-law (and possibly hard) intrinsic 
TeV spectrum over ~2 decades is generally needed (e.g. 
1ES 1426+4280). 

Figure 2 shows the effect of a change of the EBL 
normalization "P" and of source redshift: Ar (the differ- 
ence between observed and intrinsic slope) increases with 
both normalization and redshift. The two dependencies 
combine, so that at higher rcdshifts an equal change of 
EBL density yields a much stronger effect (Fig. 2, right 
panel) . Redshift thus provides leverage: more distant ob- 
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Fig. 3 Effects of changing the EBL spectral shape (increasing the flux selectively at 1-2 fim or at 0.2-0.3 jj,m, left panel) 
on the attenuation factor, for fixed redshift. Center panel: an increase only at 1-2 fim strongly steepens the spectrum, since 
it increases the optical depth at ~ 1 — 2 TeV without affecting the 0.2-0.3 TeV range (since below threshold). Instead, a 
higher EBL flux in the UV range alone increases absorption at 0.2 TeV more than at 1-2 TeV, thus reducing the energy 
dependence of the optical depth (right panel). As a result, the spectrum suffers less steepening, at the price of an overall 
higher attenuation. 



jects represent therefore a more sensitive diagnostic tool 
for EBL studies, but with the limitations given by gen- 
erally lower fluxes (and thus statistics). 

Figure 3 shows instead the effects of a change in the 
EBL spectrum with respect to our template, namely in- 
creasing the flux in the 1-3 /jm range (as would be re- 
quired to match some direct estimates, and expected in 
the Pop III stars scenario) or in the 0.1-0.3 /im range. 
In the first case the steepening gets stronger, while in 
the second case it is reduced, because of a lower con- 
trast in optical depth between 0.2 and 1 TeV photons 
(though at larger absolute values). Note that a decrease 
of the UV-opt fluxes would not reduce the steepening, 
since it increases the contrast between the two energies 
while reducing the absolute values. 

It is also worth to highlight an important consequence 
of EBL absorption on blazar interpretation: if the EBL 
is high, sources at z > 0.1 would be so heavily absorbed 
that their intrinsic Inverse Compton (IC) luminosity and 
peak energy become always much higher than the ob- 
served ones. The absorption-corrected values can eas- 
ily reach IC peak energies > 3-10 TeV and luminosi- 
ties in excess of 10 47-48 ergcm~ 2 s -1 , with a Compton 
dominance of 10-100. Therefore, the bulk of their lumi- 
nosity would be emitted at high energies, alike FSRQ. 
In other words, they would represent the famous "high- 
energy peaked, high-luminosity" blazars that would in- 
validate the blazar sequence scenario [T2l[TB] (though it 
would remain the issue of why their radio luminosity is so 
low compared to FSRQ). These objects would not have 
been recognized as such counter-examples simply be- 
cause their luminosity would be largely underestimated, 
being dispersed in the intergalactic space due to 7 — 7 
absorption. Though now this seems not the case (see 
Sect. 3), this possibility still remains also with a low EBL 
level, considering objects at significantly larger redshift 
(but which are likely beyond the horizon of the present 
generation of CT). 
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Fig. 4 The HESS spectra of 1ES 1101-232, corrected for ab- 
sorption with three different EBL SEDs, as labelled (from [5], 
see ). Red: observed data. Blue: absorption-corrected data. 
The lines show the best fit power-laws to the reconstructed 
spectrum, and the corresponding shape after absorption 



3 EBL upper limits from H.E.S.S. spectra 

The intrinsic TeV spectra were reconstructed directly 
from the observed ones, using the optical depth calcu- 
lated for the assumed EBL SED. In this way no a pri- 
ori assumption on the blazar spectrum is required. The 
reconstructed spectra are compatible with a power-law 
model, but high EBL fluxes yield extremely hard pho- 
ton indices (r int < 0), for both 1ES 1101-232 (Fig. 4, and 
whole SED in Fig. 5) and H 2356-309 (see [5]). Such hard 
spectra have never been seen in the closest, less absorbed 
objects (e.g. Mkn 421 and Mkn 501, r int w 1.5 - 2.8 for 
the same EBL SEDs) and are difficult to explain within 
the standard leptonic or hadronic scenarios [4j for blazar 
emission. Assuming instead more typical values for the 
intrinsic spectrum as expected from the currently known 
blazar physics and phenomenology, namely that the true 
average VHE spectrum during the H.E.S.S. observations 
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Fig. 5 SED of 1ES 1101-232^5], for an EBL spectrum equal 
to Pl.O. At lower energies, XMM data (blue) and RXTE data 
(green) have been taken in the same epoch of the H.E.S.S. 
observations (Aharonian et a l. 2006, in prep.). Magenta: his- 
torical SAX observations |25] . 

was 7] n t > 1.5 (as obtained under most circumstances if 
the index for the particle spectrum in shock acceleration 
models is s >1.5, [H]), the EBL SED has to be scaled 
down to the level of P0.45, very close to the galaxy counts 
limit (details in [5]). This represents the upper limit for 
the EBL in the Opt-NIR range which does not require 
dramatic changes for the blazar physics or in the sce- 
narios for the underlying particle spectra. Accounting 
for galaxy evolution and for the statistical and system- 
atic uncertainties on the H.E.S.S. spectral measurement, 
the limit corresponds to < (14 ± 4) nW m~ 2 sr _1 (i.e. 
< P0.55±0.15; Fig. 6). Of course, a different limit can be 
derived according to different hypotheses, and depending 
on their overall plausibility. But only strong spectral dif- 
ferences would qualitatively change this result (see Sect. 
4): for example, even the case of a particle spectrum with 
s = 1, which has been argued as asymptotically possi- 
ble in finit-Mach-number shocks under certain conditions 
(see e.g. [T5II16] ). would raise the limit only slightly (for 
1ES 1101-232, the change of scaling factor P is approxi- 
mately AP ~ 0.34Z\r). 



4 Alternative scenarios 

Another way to avoid such hard spectra is to reduce 
the energy dependence of the optical depth by increas- 
ing absorption at low 7-ray energies, i.e. with high UV- 
Opt fluxes. However, the huge fluxes required to have 
r = 1.5 (> 300 nW m- 2 sr" 1 , see [5] and Fig. 3 left) 
are in contrast with measurements[17], even when in- 
terpreted as upper limits, and could hardly be accomo- 
dated within any reasonable cosmological model. Un- 
less considering even more exotic scenarios like violation 
of Lorentz invariance or the non-cosmological origin of 
blazars' redshift, at present the most viable alternative 
is that such hard spectra are a real, newly discovered fea- 
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Fig. 6 Upper limit for the EBL (black marked region), cor- 
responding to the assumption of an intrinsic 1ES 1101-232 
spectrum with photon index F > 1.5. It is obtained from 
the P0.45 curve, when accounting for galaxy evolution effects 
(details in [5j). 

ture of the blazar emission. Possible mechanisms have al- 
ready been envisaged[4]. Bulk-motion Comptonization in 
the deep Klein-Nishina regime of a narrow-band photon 
field (such as a Planck-type distribution) by cold plasma 
with a very large Lorentz factor may lead to very hard 
spectra with a sharp pile-up, reproducing spectra like 
the ones in Fig. 4,5 (see [!])• A "pile-up" (i.e. relativis- 
tic Maxwellian) particle distribution seems also expected 
as outcome of turbulent acceleration (see e.g. [TSUll?] 
and refs therein). In this case, however, even suppressing 
cooling and considering a purely monoenergetic distribu- 
tion (or a distribution with a similarly sharp low energy 
cut-off at 7mm > 10 5 [ID]), the resulting synchrotron and 
SSC spectrum would be at most cx j/ 1 / 3 , i.e. r > 0.66 
( [2D]). If so, the EBL P0.45 upper limit would be shifted 
to ^P0.72. Note that to make this scenario feasible it is 
necessary to always hide the hard synchrotron emission 
of these particles below a second, broader component 
which dominates the blazar emission at low energies, in 
order to account for the simultaneous SED (Fig. 5). 

However, even if theoretically feasible, also this al- 
ternative seems highly unlikely. If intrinsic, such pile-up 
features should be directly visible in the observed VHE 
spectra of closer, less absorbed HBLs like Mkn 421 and 
Mkn 501 (if 7] n t was « 0, they should show r o b s < 0.5). 
This is in contrast with observations, unless assuming a 
dependence of the source parameters on redshift such 
that the corresponding features always disappear due 
to EBL absorption. It is difficult to justify such a fine- 
tuning on a relatively small redshift range, even assuming 
an evolution of blazar properties. Though the sample of 
objects is still too limited to settle this issue definitively, 
growing data seem to corroborate this conclusion. Figure 
7 shows the observed photon indices of all TeV blazars 
detected so far as a function of redshift, in comparison 
with the expected index for a r = 1.5 intrinsic spectrum. 
If the EBL density is high, the very strong dependence 
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Fig. 7 Observed hardest photon indices of the VHE spec- 
tra of all detected TeV blazars as a function of redshift. 
The lines give the expected observed (i.e. absorbed) pho- 
ton index of a source with an intrinsic spectrum t = 
1.5, when fitted in the range 0.2-1 TeV, for different EBL 
SEDs (as labelled). The dashed line corresponds to an 
EBL SED as low as the galaxy counts limit, and including 
galaxy evolution effects (according to [11]). The distance of 
the points over each curve represents visually the intrinsic 
hardness of the source spectra. Red: H.E.S.S. data. Blue: 
MAGIC data. Black: past results. The sources are (from 
left to right): Mkn421, Mkn501, 1ES2344+514, Mkn 180, 
1ES 1959+650, PKS 2005-489, PKS 2155-304, 1ES 1426+428, 
H 2356-309, 1ES 1218+304, 1ES 1101-232, PG 1553+113 (de- 
tails in 16;, Table 1). 

on redshift would require a dramatic change of properties 
between TeV blazars at z ~ 0.2 and z < 0.1, with all 4 
farther objects showing such hard features (which would 
thus seem rather common) and none of the 8 closer ob- 
jects. There seems to be no reason why the jet emission 
mechanisms in blazars should know the level of the EBL 
with such precision. A low EBL level instead, close to 
the galaxy counts limit, does not impose such difference, 
making the TeV properties similar among different HBLs 
in this redshift range, as observed in all the other energy 
bands of the SEDs. 



5 Observational tests 

A straightforward step to test this conclusion is to look 
for independent evidence that HBLs can indeed have 
much harder TeV spectra than assumed. This can be 
done focussing in particular on lower (z < 0.1) and 
higher (z > 0.3) redshift objects. At low redshift ab- 
sorption is less, so hard features can be directly mea- 
sured in the observed TeV spectra (e.g. r o b s < 1). New 
detections and a monitoring program will increase the 
chances to catch an object in these special states, should 
they occur seldom. Observations of high-redshift targets, 
in addition to extending the coverage of the r — z plane 
(Fig. 7), could reveal objects whose reconstructed spec- 
trum is still very hard even with an EBL flux as low 
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Fig. 8 Left, full line: Primack 2001 model[TT|, used to cal- 
culate the absorption properly including galaxy evolution ef- 
fects. The original model was scaled down 0.8 x to match the 
new EBL upper limit (dashed line). Right: corresponding at- 
tenuation factor for a source at z=0.5 (black), z=0.116 (red) 
and z=0.186 (blue). 

Table 1 High-z TeV candidates for MAGIC/ VERITAS 



Name 


z 


Sample 


1ES 1106+244 


0.460 


Einstein Slew 


RX J1326. 2+2933 


0.431 


RASS 


1ES 0219+428 


0.444 


Einstein Slew 



as the galaxy counts limit. This would provide evidence 
for intrinsically hard spectra, though in such case the 
possibility of exotic scenarios invoking a modification of 
the physics of the 7-7 interaction would not be formally 
excluded. With a low EBL level, objects at z=0.4-0.5 
should be well within the possibilities of the current CT 
generation. Figure 8 shows the optical depth expected 
for a source at z=0.5 using the EBL model by [TT], when 
scaled down to match the EBL upper limits. This model 
was adopted to properly include galaxy evolution effects, 
which become important at larger redshifts. Though the 
spectrum is expected very steep, up to 300 GeV the at- 
tenuation is not dramatically higher than for the already 
detected TeV sources, and only a factor ^2 around 150 
GeV. Sources at such distances therefore could be de- 
tectable by current CT also outside exceptional (and 
rare) high-flux states. 

The best targets for such observations are HBLs with 
the appropriate declination to be observable by CT at 
zenith (to take advantage of the lowest possible energy 
threshold) , and with high X-ray brightness (indicating a 
large number of TeV electrons). Using the BLLac sam- 
ples and criteria described in [2T], a group of 3 objects 
best matching these requirements have been found in 
the northern hemisphere (Fig. 9), suited for MAGIC and 
VERITAS. They are therefore suggested for observations 
(Table 1). A similar program for the southern hemisphere 
has been proposed for H.E.S.S. . 



6 Conclusions 

Showing an unexpected hardness for their redshift, the 
combined H.E.S.S. spectra of the HBLs 1ES 1101-232 
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Fig. 9 BL Lac objects in t he radio (5 GHz) and X-ray (1 
keV) vF(y) plane (details in [21]). Sources belonging to differ- 
ent samples have different symbols, as labeled. Black crosses: 
new TeV detections since the proposal of the X-ray/radio 
Flux selection criterium (dotted rectangle). Black squares: lo- 
cation of the proposed high-z sources for MAGIC/ VERITAS 
(Table 1). If at redshift ~ 0.1, their location would be in the 
upper part of the rectangle (a change of redshift moves an 
object approximately along the lines of constant qrx). 



and H 2356-309 (now supported by two further TeV sources 
provide strong circumstantial evidence for a low density 
of the EBL. With minimal assumptions on the blazar 
properties, the H.E.S.S. spectra allow the derivation of 
an upper limit on the EBL in the Opt-NIR range which 
is very close to the lower limit given by the integrated 
light of resolved galaxies. This result indicates that the 
EBL at these wavelengths is strongly dominated by the 
direct starlight from galaxies, and thus excludes a large 
contribution from other sources, in particular Pop III 
stars. At the same time, this result strongly reduces the 
uncertainty on the reconstruction of the blazar intrinsic 
spectra (up to few TeV), thus allowing a more accurate 
modelling and understanding of the overall SED. How- 
ever, additional detections and monitoring by the new- 
generation CTs are needed to better sample the blazars' 
behaviour at VHE, to further test this conclusion and 
possibly to pin down the diffuse background also above 
a few microns. 
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